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Abstract 

The classical nucleation theory of Becker, Doring and Zeldovich is compared with 
the Langer coarse-grained field approach to the nucleation phenomenon. Both for- 
malisms have been applied to the condensation from a supersaturated vapor. It is 
shown that the nucleation rate derived in the classical theory can be expressed in a 
form equivalent to that of the field nucleation theory. This equivalence serves as an 
explanation of the puzzling fact that the numerical predictions of both theories for 
condensation of Xe and CO2 are almost identical though the standard analytical 
expressions for the nucleation rates are different. The results obtained can help to 
link the theories of nucleation and their approximations. 
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First order phase transitions, which include the thoroughly studied condensa- 
tion of a vapor, boiling of a liquid or crystallization of a melt, as well as the 
hypothetical transition from a hadronic phase to the quark— gluon plasma, 
play an important role in modern science and technology. The theory of first 
order phase transitions is based on the homogeneous nucleation theory which 
aims to determine the rate of the decay of a metastable state, i.e. the rate 
of formation of nucleating clusters within the initially homogeneous nonequi- 
librium state. This is the first important assumption of the nucleation theory 
which implies that the conversion of a metastable state into a thermodynam- 
ically stable phase should proceed via the formation of clusters of the new 
phase. As a typical example one may consider the condensation of small drops 
of liquid from the vapor rapidly quenched into the metastable region. 

Clusters of various sizes are produced spontaneously due to the fluctuations in 
a homogeneous medium. Clusters with sizes larger than a certain size, called 
critical, will grow into a new phase, and thus the phase transition process de- 
velops. A description of this process requires a kinetic analysis of the evolution 
of the clusters. The phenomenological nucleation theory, often referred to as 
the classical or Becker— Doring—Zeldovich (BDZ) theory, was formulated in 
the pioneering works by Becker and Doring [1], and Zeldovich [2]. It is ex- 
plained in detail in many textbooks and reviews, see e.g. Refs. [3-9]. Despite 
the phenomenological origin, the classical theory appears to predict nucleation 
phenomena very well. 

A fleld theoretical approach to the problem of nucleation has been proposed 
by Langer [10-12] who extended the generalization of the BDZ theory for a 
system of arbitrary many degrees of freedom [13] to the field theories. This 
formalism has been applied to several nucleation phenomena, including the 
liquid— vapor phase transition [12]. In spite of all the statistical and hydrody- 
namic corrections, it appears that the nucleation rate calculated within the 
coarse-grained field theory differs only slightly from that of the BDZ the- 
ory, while the expressions for the rates are different. This surprisingly good 
agreement between the two theories leads to the assumption that, although 
the coarse-grained field theory is a more accurate and advanced treatment of 
the nucleation phenomena than the classical theory, the latter compensates 
somehow correctly for all degrees of freedom [7,8,12]. 

In the present paper we would like to show that the nucleation rate derived 
in the BDZ theory may be obtained under certain assumptions in the form 
which is identical to that of the Langer theory. We will start from a brief review 
of the classical nucleation theory. Then, after discussion of the formalism of 
the coarse-grained field theory, the results of both models will be applied to 
condensation from a supersaturated vapor. The comparison will be done in 
terms of the Langer theory. 
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The metastable state typically is assumed to be a mixture of molecules and 
clusters of the new phase containing n molecules. Let Nn{t) be the average 
number of clusters of size n at time t. Effects of the direct interactions between 
the clusters are disregarded. Therefore, the clusters may change their size only 
due to an evaporation— condensation mechanism, in which a cluster of size n 
grows or shrinks by condensation or loss of a single molecule. This is the basic 
assumption of the classical nucleation theory. The rate of change of the number 
of clusters of a given size is determined by the master equation 

— = Jn - Jn-1 , n>2 , (1) 

where 

J„ = y+_i7V„_i(i)-y-7V„(i) (2) 

is a current in a size space, i.e. it is the rate at which clusters of size n — 1 grow 
to clusters of size n. Here g^-i is the rate of the molecule condensation to a 
cluster of size n, and 9,^7 is the rate of evaporation. According to the theory 
of thermodynamic fluctuations the equilibrium distribution of clusters obeys 
the Boltzmann (or Gibbs) distribution 

K'-NS'..,{-^) . (3) 

where AG„ is the minimum work needed to form a cluster of size n, is 
Boltzmann's constant, T is the temperature, and A^q* is a pre-exponential 
factor. A maximum in AG„ is reached at dAGn/dn = corresponding to 
a cluster of critical size, Uc- Clusters smaller than the critical cluster tend to 
shrink, clusters of the critical size are in a metastable equilibrium, and clusters 
larger than the critical cluster tend to increase their volume. 

The equilibrium cluster distribution function N^'^ is a time independent solu- 
tion of the master equation (1). The current J„ is zero for this distribution, 
therefore, one may write the equation of detailed balance 

As shown by Zeldovich [2] , under the treatment of n as a continuous variable 
the master equation (1) may be approximated by the Fokker— Planck kinetic 
equation 



dt dn dn 



kBT dn ' ' ^" dn 



(5) 
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or, in terms of the ratio Nn{t)/N^^, 



dNn{t) 
dt 



d_ 
dn 



dn I 



(6) 



The classical theory is based on a nonequilibrium, steady-state solution of 
Eq. (5), corresponding to a continuous phase transition. Then Jn = s — 
constant, which is called the nucleation rate and measures the number of 
clusters passing through the critical range per unit time per unit volume. 

The boundary conditions for small and large n are chosen in the BDZ theory 
as follows, 

= l , ltoM)=o. (7) 



A number of small clusters is assumed to be close to its equilibrium value due 
to the rapid increase of the fluctuation probability with decrease of cluster size. 
It is assumed also that clusters larger than that of critical size are removed 
from the system. The steady-state solution of Eqs. (6) and (7) reads 

Ml = . 7 r8^ 

I 9'nKr 



where 

dn' 



n' 



(9) 



Since the function N^''- = Nq'^ exp^—AGn/ksT) has a very sharp maximum 
at ric, the integral in (9) can be evaluated by the saddle-point method, i.e. by 
expanding of AGn around nc and extending the integration with respect to 
n — nc from — oo to oo. With these approximations, the result for the steady- 
state nucleation rate may be written in the general form 

s = g^NllZ , (10) 



containing the Zeldovich factor 
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Here the number of clusters that have passed through the critical range is 
expressed via the equilibrium number of clusters of critical size. This is the 
principal result of the classical theory. 

The accretion rate of the molecule condensation to the surface of the critical 
cluster is derived from the kinetic theory of gases to be 

_ pSc 



containing the gas pressure p, the surface of the critical droplet Sc, and the 
average mass of a single molecule mi. On the other hand may be obtained 
from the solution of rather general equations describing the diffusion growth 
of dense clusters. We will come to this point later during the discussion of 
vapor condensation. 

An extension of the classical theory to a system of arbitrary many degrees of 
freedom was worked out by Landauer and Swanson [13]. Their method was 
applied by Langcr to the field nuclcation theory [11]. The BDZ and the Langer 
theories postulate that unstable fluctuations of thermal origin, e.g. clusters of 
critical size, lead to the decay of the metastable state. However, in the field 
nucleation theory the critical cluster needs not necessarily be a physical object 
[14], but rather characterizes a certain saddle-point configuration. In contrast 
to the classical theory, only the critical cluster and its shape deformations are 
required for the field theory of nucleation. The advantages and shortcomings 
of both models are discussed in detail in [7] . 

The Langer theory starts with the consideration of a classical system with A/" 
degrees of freedom (coordinates and momenta) rji, i = 1, . . . ,J\f. The probabil- 
ity distribution functional p{{r)i},t) is assumed to obey a continuity equation 
in {77}-space, namely 



which is analogous to the multi-dimensional Fokker— Planck equation. The 
probability current Jj is given by 

^. = -gM«gp+..T^). (14) 



containing a generalized mobility matrix M^j, and a coarse-grained free energy 
functional F{rj}. For the case of Jj = the equihbrium solution of Eq. (13) 
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reads 



PeqiVi} « exp 



( 



(15) 



The basic idea of the field nucleation theory is that the states of metastable 
and stable equilibrium lie in the vicinity of the local minima of F{ri,i} in the 
{?7}-space, which maximize Peq{fli}- The phase transition occurs when the 
configuration {77^} moves from the vicinity of a metastable minimum to the 
vicinity of a stable minimum. To develop into the stable state the configuration 
must overcome the potential barrier. It is most likely for the trajectory of the 
system to pass across a small area around the intermediate saddle point of 
the functional F, say {r^}. This saddle-point configuration {f^} corresponds 
to the critical cluster of a condensing phase in the BDZ theory, although it 
may not be a physical object. The rate of the decay of the metastable state is 
determined by the steady-state solution of Eqs. (13) and (14), corresponding 
to a finite probability current across the saddle point from the metastable to 
the stable minimum of F. 

To obtain this solution Langer [11] introduced a set of new variables ^k: k — 
1, . . . ,J\f, whose fluctuations are slow compared to the characteristic time of 
the evaporation— condensation process. These variables should contain the or- 
der parameter. They are related to the r/'s by virtue of the orthogonal trans- 
formation D: 



where, near the saddle point {^} = {0}, 

m = ^+^EA,e|+--- (17) 

Here F = F{f]}, and Xj are the eigenvalues of the matrix Mij = {d'^ F / driidirij) ,^=f^. 
By the definition of the saddle point at least one of the eigenvalues \j must 
be negative. Following Langer we denote this eigenvalue as Ai. Then in terms 
terms of the ratio cr({{}) = p({C})/Peg({C}) ^^e requirement of a steady state 
current over the potential barrier has the form 



which is a straightforward generalization of Eq. (6). The stationary probability 
distribution functional p({C}) obeys the boundary conditions that, near the 



(16) 




(18) 
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mctastablc state, p{{C}) is very close to the equilibrium distribution PeqiVi}^ 
and beyond the saddle point should vanish rapidly. The solution of 

Eq. (18) which satisfies these boundary conditions becomes 



where m is a linear combination of the ^^s 

u^Yl UiCi . (20) 

i 



The coefficients Ui are solutions of the eigenvalue equation 

A, J2 MmU, = -kU, , (21) 

i 



which describe the growth rate k of the single unstable mode at the saddle 
point. To obtain the desired nucleation rate one has to integrate the fiux 
Ji over any plane containing the saddle point and not parallel to After 
the cumbersome but straightforward calculations the steady state solution of 
Eq. (18) may be written in the form [11] 



where 



Here AFg = F{fi) — F{rio) is the excess free energy required to form a cluster of 
the critical size. The dynamical pref actor k describes the cluster growth at the 
saddle point, the statistical prcfactor Qq is often called the generalization of 
the Zeldovich factor, and V denotes the volume of the saddle point subspace. 
Although Eq. (22) has a structure similar to Eq. (10), it is asserted that the 
pref actors in both equations should be different [12]. However, up to now there 
is no crucial experiment which can distinguish the predictions of Langer theory 
from those of the BDZ theory. 

To compare the models let us consider the condensation of a supersaturated 
vapor. The minimum work needed to form a spherical droplet of radius it! in 
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the thin wall approximation is simply 

AF{R) = -—R^Ap + AnR^a , (24) 



where Ap is the difference between the pressures inside and outside the droplet 
and a is the surface tension. The critical radius is given by the Laplace formula 
Rc = 2(7/Ap, and near this value AF{R) can be represented by harmonic 
approximation 

AF{R) = ^naRl - ina{R - R^f . (25) 



It is convenient to rewrite the excess of the free energy in terms of the reduced 
radius r — R/Rc and the ne 
by the authors [15] recently, 



/ 47r(j 

radius r = R/Rc and the new similarity number \z — Rc\ — introduced 



kf^T 



AF{R) = - kBT{r - ifXl (26) 



Prom the definition (17) it follows that ^ = f — the free energy of the 
critically large droplet is F = ksTX^/S and the only negative eigenvalue Ai is 
Ai = 2kBT\\. Therefore, the nucleation rate predicted by the Langer theory 
becomes 



The statistical prefactor is obviously determined by the formula 



Substituting Eq. (26) in Eq. (11) we get for the Zeldovich factor 



/ \R, 



and 

^ ^ ^ (30) 
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The factor Z has the physical meaning of a measure of the width A of the 
Gaussian distribution by which the actual distribution p^^[K) is approxi- 
mated, namely Z = A"^. Therefore, the statistical prefactor Qq, which mea- 
sures the volume of the saddle point in the phase space, can not be considered 
as a simple generalization of the Zeldovich factor. Instead, it is related to the 
Z factor by 

Qo = ^ . (31) 



The nucleation rate derived in the BDZ theory is 



where one has to evaluate the rate g{Rc)- Prom the significance of the current 

g d^F 

J, the drift coefficient — acts as a velocity in size space [6], dR/dt, 

kbT oR 

which can be determined independently by the solution of the macroscopic 
equations for the growth of a droplet. We will consider the diffusion model 
of droplet growth which assumes that the droplet of a new phase grows due 
to the diffusion flux of matter through its interface. Then, the problem of the 
liquid— gas phase transition is similar to the case of precipitation of a sub- 
stance from a supersaturated solution or separation in binary fluids. The only 
difference is that one has to choose the order parameter as the local entropy 
density in a liquid— gas system, and as the local concentration in a weekly 
supersaturated sohition [17]. For the latter case the spherically symmetric 
concentration distribution c(r) aroimd a droplet of radius R is given by the 
steady state solution of the diffusion equation 

DV'cir) = ^ = , (33) 



where D is the solute diffusion coefficient. We obtain 

R 

c(r) = c- (c-cor)— , (34) 



where 



R 

cqr = Co H (c - Co) (35) 

r 



is the equilibrium concentration of the solution at the surface of a droplet of 
radius R, Cq is the equilibrium concentration of the solution above the planar 
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surface, and c is the given value of the concentration of the supersaturated 
solution. Substituting Eq. (35) in Eq. (34) and identifying be as the supersat- 
uration c — Cq in the system, we get 



(36) 



dc 

Since the diffusive flux into the droplet i(r) = £)— at a droplet surface is 

or 

equal to the growth rate of the droplet, one may write 



dR 
'dt 



DSc 



^ Rc 
R 



(37) 



Linearizing this equation around the critical radius and taking into account 
that, by definition of the dynamical prefactor, R — Rc — exp {—Kt), we have 
finally 



dR , , /RA^ „ , 



(38) 



where \k\ 



DSc 



. This radial dependence of the droplet growth rate is valid for 

the liquid— gas phase transition also. Similar result was obtained in Ref. [16] 
by virtue of the hydrodynamic Kotchine equations describing the enlargement 
of a droplet. Thus, for the critically large droplet we have 



9AF 
OR 



Rc 



\K,\ 

2A| 



(39) 



Substitution of Eq. (39) in Eq. (32) yields 

jBoz = ^ , (40) 

leading again to the result (27). 

For the sake of simplicity we considered the process where the one order pa- 
rameter was sufficient for the description of nucleation. Obviously, our analysis 
may be extended to the system of arbitrary many degrees of freedom. We pro- 
posed a direct method allowing one to link the classical theory to the modern 
field theory of nucleation. It is very important for the comparison that the 
underlying dynamics of the process, that governs the cluster growth, should 
be the same. Trying to compare both theories Langer and Turski found [12], 
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that the numerical predictions of the critical supersaturation in Xe and CO2 
given by the field micleation theory are very close to those of the classical 
theory. But the rate of molecule accretion on the surface of the critical clus- 
ter, g{Rc), was calculated from simple geometric considerations of the kinetic 
theory of gases. Therefore, the analytical expressions of the nucleation rates 
were different, and it was difficult to find conformity between them. 

Recently [18] the Langer theory of nucleation has been applied to calculate 
quark— gluon plasma formation in heavy ion collisions at the energy of the 
Brookhaven accelerator. Then the nucleation rate derived in the field theory 
has been compared with that of the classical theory. Again it turns out that 
both prefactors are about the same order of magnitude. 

The classical nucleation theory is attractive not only because it predicts nu- 
cleation rates in terms of macroscopic quantities. It can be easily generalized 
to time-dependent nucleation phenomena. Also, the classical theory can be 
extended to the temperature and curvature dependences of the interfacial en- 
ergy, to the deviations of the shape of the clusters from the ideal spherical 
shape, etc. In spite of these advantages, and although the classical theory ap- 
pears to be consistent with the experimental data, it is implied explicitly and 
implicitly that one must consider the classical theory only as a crude approx- 
imation. This problem becomes clearer now. Using the macroscopic diffusion 
model for the description of a spherically growing droplet of liquid we have 
shown that the classical nucleation rate can be derived in a form identical to 
that of the field nucleation theory. The results obtained in this paper support 
the validity of the phenomenological classical theory on the basis of the more 
accurate Langer first-principles theory of nucleation. A further comparison 
between these theories may be useful for the development of a more rigorous 
nucleation theory. 
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